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We discuss a constant contribution to meson correlators at finite temperature. In the deconfine- 
ment phase of QCD, a colored single quark state is allowed as a finite energy state, which yields 
to a contribution of wraparound quark propagation to temporal meson correlators. We investigate 
the effects in the free quark case and quenched QCD at finite temperature. The "scattering" con- 
tribution causes a constant mode in meson correlators with zero spatial momentum and degenerate 
quark masses, which can dominate the correlators in the region of large imaginary times. In the free 
spectral function, the contribution yields a term proportional to u)5(uj). Therefore this contribution 
is related to transport phenomena in the quark gluon plasma. It is possible to distinguish the con- 
stant contribution from the other part using several analysis methods proposed in this paper. As 
a result of the analyses, we find that drastic changes in charmonium correlators for \ c states just 
above the deconfinement transition are due to the constant contribution. The other differences in 
the Xc states are small. It may indicate the survival of \ c states after the deconfinement transition 
until, at least, 1.4T C . 

PACS numbers: 12.38.Gc,12.38.Mh 



I. INTRODUCTION 

In the imaginary time formalism of finite temperature 
field theory, temporal correlators are related to several 
dynamical properties through a prescription of analyti- 
cal continuation. Especially the imaginary time corre- 
lator can be described by the same spectral function as 
that of retarded and advanced green functions. There- 
fore the spectral functions extracted from the imaginary 
time correlators, which are able to be calculated in lat- 
tice QCD, enable us to investigate important physics at 
finite temperature, e.g. dilepton rates [lj, light mesons 
0, heavy quarkonia [j, 3, [E IE Bi 3 , glueballs @, trans- 
port coefficients [Tol. [ill Tl3j . 

One of the most important quantities at finite temper- 
ature QCD is the spectral function of heavy quarkonium, 
which plays the key role for understanding the quark 
gluon plasma (QGP) formation in heavy ion collision ex- 
periments e.g. the RHIC experiment at Brookhaven Na- 
tional Laboratory. Recent studies on the spectral func- 
tion of charmonium above T c suggest that hadronic exci- 
tations corresponding to J/ijj may survive in the decon- 
finement phase till relatively high temperature [J, [H, H, 0, 
8]. Such results of strongly interacting QGP may affect 
the scenario of J / ip suppression [ID, [3] ■ Therefore the 
determination of the accurate dissociation temperature 
of J/ifi is required for many phenomenological studies. 

As somewhat different from the "direct" J/ifj suppres- 
sion, recently an "indirect" J/ip suppression was pro- 
posed flol fig . The total yield of J/ -0 is not only from 
the direct production, about 40% of the J/tp's that come 
from the higher states ip' and \c \U\ ■ If the higher states 
dissociate at lower temperatures, a part of the J/ip sup- 
pression may be observed in experiments at a lower tem- 
perature than that of J/tp. Therefore the study of ip' 
and Xc states and their dissociation temperatures is a 
very interesting subject for the physics of QGP. Several 
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FIG. 1: A sketch of quark line diagrams for a meson-like 
correlator in a system with finite temporal extent. Vertical 
lines show the boundaries in temporal direction. 



groups have already investigated the Xc states, and their 
results indicate that the lowest peak of the spectral func- 
tion, corresponding to the Xc state, may dissociate just 
after the transition temperature 0, 0, H} ■ On the other 
hand there is no lattice QCD study on the ip' states above 
T c , because it is difficult to extract information of non- 
lowest states such as ip' from the correlator on the lattice 
at finite temperature. 

In such studies of temporal correlators at finite tem- 
perature, one has to take account of a characteristic con- 
tribution caused by the finite temporal extent. Usually a 
meson correlator is interpreted by a diagram with quark 
and anti-quark propagators like it is sketched in Fig.QJa). 
However in the case of a system with a finite temporal 
extent, a wraparound (scattering) contribution, as shown 
in Fig. (Hb), has to be included as well. It is similar 
to the situation at zero temperature e.g. the study of 
two pion correlators [18| , and pentaquark correlator [l9( . 
Because these correlators include multi-hadron states as 
intermediate states, and each state is also allowed as a 
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color singlet state, which can have a finite energy state in 
the confined phase. At zero temperature the calculation 
should be performed, in principle, with infinite temporal 
extent, therefore, the wraparound effect is not physical. 
The contribution is usually removed using e.g. Dirichlet 
boundary conditions or appropriate analyses jl8l |. 

In finite temperature calculations using the imaginary 
time formalism, the temporal extent is determined by 
(the inverse of) the temperature, and the quark fields 
have to have anti-periodic boundary conditions. In the 
confinement phase, an expectation value of the diagram 
in FigQJb) vanishes because it is related to a single quark 
propagation, which is exponentially suppressed with its 
infinite energy due to the confinement. On the other 
hand, in the deconfinement phase, the contribution may 
yield a finite expectation value. This is similar to the dis- 
cussion of the Polyakov loop expectation value. For ex- 
ample in the free quark system, when one adopt a meson- 
like operator with degenerate quark masses and zero spa- 
tial momentum, the correlator of the diagram in FigQJa) 
behaves like exp(— m q t) x exp(— m q t) = exp(—2m q t), 
where m q is the quark mass (or the energy of a single 
quark state). (Here we note that equations are writ- 
ten with dimension-less variables throughout this paper, 
if there is no comment on the units.) On the other 
hand, the diagram in FigfTJb) behaves like exp(— m q t) x 
exp(— m q {L t — t)) = exp(— m q L t ), where L t is the in- 
verse of the temperature. The latter diagram provides a 
constant contribution to the correlator, and the contri- 
bution may dominate in the large imaginary time region 
as a zero energy mode. Due to the contribution, cor- 
relators of such meson-like operators, e.g. charmonia, 
may be drastically changed just after the deconfinement 
transition. The latter type diagram, in general, yields 
t-dependent contributions which depends on a difference 
between single quark energies for each quark propaga- 
tor. Even in non-degenerate quark masses or non-zero 
spatial momentum cases, the diagram can yield the con- 
stant contribution when the energies are identical. How- 
ever we do not discuss the case in this paper. The finite 
momentum case has been discussed in Ref. [2(| at high 
temperature limit. 

This paper is organized as follows. In Sect. [H] we dis- 
cuss the correlators of meson-like operators and its spec- 
tral functions in a system of free quarks with finite tem- 
poral extent. Some numerical results are also presented. 
In this section we also discuss the physical interpretation 
of the constant contribution and some analysis methods 
to avoid this contribution. In Sect. IIIII we show numeri- 
cal simulations in quenched QCD on anisotropic lattices, 
and demonstrate effects of the constant contribution in 
the interacting case and the analyses to avoid the contri- 
bution. As a result of these analyses we show that the \c 
state yields a small change except for the constant mode 
in its spectral function after the deconfinement transition 
like the J/tp state at least up to T = IAT C . At last we 
discuss results on the constant contribution. 



II. THE CONSTANT CONTRIBUTION IN THE 
FREE THEORY 

A. Meson-like correlators in the free quark case 

In order to investigate the constant contribution as 
mentioned in the Introduction, we first consider the free 
quark case of QCD, in which the constant contribution 
can be easily calculated. Here we define the meson- 
like correlators with quark bilinear operators, Or (if, t) = 
q(x,t)Tq(x,t), in the free quark case this gives for the 
correlator, 



C(t) = J2(Or(x,t)Ol(0,0)), 



(1) 



where T are appropriate 4x4 matrices, i.e. 75, ji, 1, and 
7^75 for pseudoscalar(Ps), vector(V), scalar(S), and ax- 
ialvector(Av) channels respectively. In this section and 
later we always adopt the anti-periodic boundary condi- 
tion for quark fields in temporal direction to suppose a fi- 
nite temperature case, and periodic boundary conditions 
in spatial directions. In this paper we calculate correla- 
tors for degenerate quark masses, constructed from bilin- 
ear operators, with vanishing spatial momentum as the 
simplest case, which is aiming at studies of charmonium 
at finite temperature discussed in Sect. HTT1 The spectral 
function of the correlator is defined by 



C(t) 



dtupr(uj)K(uj, t), 



cosh {u{If - t)) 
sinh(w%-) 



(2) 



and has been calculated in Ref. [2CJ, |2l| in the high tem- 
perature limit. 

In the spectral representation, the constant contribu- 
tion provides a term proportional to ujS(lu), which corre- 
sponds to a constant mode in the meson-like correlators. 
This term is an odd function in frequency uj, which is 
one of the basic properties of bosonic spectral function. 
From the results in [2(J [HJ one can find several char- 
acteristic properties of the constant contribution. First 
there is no constant mode in the Ps channel, however the 
other correlators have the constant mode, which comes 
from the scattering contribution. The correlator in the 
S channel has also no constant mode at vanishing quark 
mass, where the chiral symmetry is restored. Secondly, 
in the V channel (and the Av channel in the chiral limit), 
there is no constant contribution coming from the zero 
relative momentum process. On the other hand, there is 
one in the Av and S channels (and it disappears in chiral 
limit) 1 . These properties are very important especially in 



1 One can find this property from the Table I of Ref. [2lll . No 
constant mode and zero relative momentum correspond to d^' = 
and d = (if quark mass is finite) respectively. 
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a finite volume calculation. When the minimum nonzero 
momentum due to the finite volume is large, the contribu- 
tion of the zero relative momentum processes dominates 
in the correlators. Furthermore, a two quark state with 
zero relative momentum is forbidden in the P-wave state. 
Therefore the relative contribution of the constant mode 
increases exponentially as the spatial volume decreases in 
the P-wave state at finite quark mass. As a result, we can 
expect a very different behavior of correlators between S- 
wave and P-wave states in the case of finite quark mass 
and finite (but relatively small) volume. Thirdly, corre- 
lators for each channel coincide at the midpoint t = L t /2 
as mentioned in Ref . [20| . These characteristic properties 
will be confirmed also by the numerical results presented 
in the next section (|IIB|) . 

It is important to consider physical interpretations of 
this contribution. The term lu5(lo) remains in the con- 
tinuum form and infinite volume results [20L [2l| . while 
it disappears at zero temperature, i.e. infinite temporal 
extent. Therefore the contribution is caused by a kind of 
physical thermal effect. The origin of the effect could be 
considered as a scattering with the thermal bath (Landau 
dumping) which is related to some transport coefficients. 
From the Kubo formula, for example, a derivative of the 
spectral function in the V channel, pv{w) which is de- 
fined by Eq. ([2]) using V = ji, is related to the electrical 
conductivity a [111 ]. 



ld_ 

6 duj 



(3) 



Let us also note on the dilepton rate. The dilepton 
rate is extracted from the spectral functions in the vector 
channel of spatial (pv{oj) with T = 7,) and temporal {pv 
with r = 70) components, 



P dilepton (u>) = p v (u) + PvM- 



(4) 



pv (w) is the charge susceptibility and it has only the 
ujS(lo) term which is equivalent to that of py(w) with the 
opposite sign. Therefore the constant contribution in the 
dilepton rate cancels at least for the free quark case. 



B. Free quark calculations on a lattice 

In this subsection we present numerical calculations of 
the meson-like correlators for free quarks on lattices. The 
calculations are performed on isotropic x N t = 16 3 x 
32 lattices, and the volume dependence is discussed using 
results for 96 3 x 32 lattice. The free quark is described 
by the Wilson quark action with a bare quark mass of 
vriq = 0.2. The boundary conditions are anti-periodic in 
temporal direction and periodic in spatial directions. 

Figure [2] shows correlators (Eq. ([I])) for each chan- 
nel. One can clearly see a constant contribution in the 
S and Av channels. The qualitative difference between 
the (Ps,V) and (S,Av) channels and the coincidence of 
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FIG. 2: Meson-like correlators for each channel. These are 
free quark calculations on a 16 x 32 lattice. The scale is 
logarithmic on the vertical axis. 
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FIG. 3: Effective masses of meson-like correlators for each 
channel. These are free quark calculations on a 16 3 x 32 lat- 
tice. 



correlators at the mid point t = N t /2 can be explained 
as mentioned in the previous section. 

Figure [3] shows the effective masses of the correlators 
in Fig. [2] The effective mass is defined by correlator at 
successive time slices. 



C(t) 



cosh [m eS (t)(f-t)] 



C(t + 1) cosh [m eff (t) (f 



t-1 



(5) 



The effective masses of the meson-like correlators with 
free quarks should approach the energy of the two quark 
state without momentum for S-wave states (Ps and V 
channels) and with a minimum momentum for P-wave 
states (S and Av channels) except for the zero energy 
mode. In the Wilson quark action, the free quark prop- 
agator satisfies the dispersion relation, 



cosh E q (p) = 1 



sm z { Pi ) + (to, +P72) 



2(1 



•p2/2) 



(6) 



where pi = 2 sin(pj/2) and E q (p) is the energy of a quark 
with momentum pi. Using Eq. ([S]), one can calculate the 
lowest energies of the two quark states in S-wave and 
P-wave states. The values obtained for N s — 16 and 
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FIG. 4: Meson-like correlators for each channel. These are 
free quark calculations on a 96 3 x 32 lattice. The scale is 
logarithmic on the vertical axis. 
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C. An analysis to avoid the constant contribution 
and Polyakov loop sectors 

Since the lattice QCD is formulated in Euclidean 
space-time, we can isolate the lowest state contribu- 
tion from correlators in the large Euclidean time region, 
while studies of non-lowest states are not so easy. When 
one want to know the information about nonzero energy 
states, the constant mode makes it difficult. The lowest 
scattering process yields only a constant mode in the cor- 
relators considered in this paper, and does not affect the 
other dynamics at higher energy, to > 2m q , at least in the 
free quark case. Even in the case of interacting QCD, if 
the width of the peak structure in the spectral function 
at ijj — is sufficiently small, e.g. u -C T, the con- 
tribution yields only a constant mode in the correlators 
[ToL |22| . Therefore, in order to investigate the nonzero 
energy states in the correlators, it is useful to remove 
the constant mode from the correlators in an appropri- 
ate analysis. In this subsection we present two analysis 
methods to avoid the constant contribution from the cor- 
relators. 

Of course the cosh + constant fit enables us to do so, 
but here we present the methods without a fit analy- 
sis. The constant mode in the correlators can be re- 
moved by the derivative of the correlator with respect to 
t [23| , where the differential correlator C (t) is defined by 
a derivative of C(t), 



C'(t) = C(t + l)-C(t). 



(7) 



FIG. 5: Effective masses of meson-like correlators for each 
channel. These are free quark calculations on a 96 3 x 32 lat- 
tice. 



Here C'(t) is equivalent to the symmetric derivative at t+ 
1/2. One can define the effective mass of the differential 
correlator, 



uiq = 0.2 are shown in Fig. [3] In the Ps and V chan- 
nels, the zero relative momentum processes have no con- 
stant contribution and their effective masses approach 
the expected values. In the S and Av channels, on the 
other hand, their effective masses approach the zero en- 
ergy level rather than the energy of two quark states. 

Next I present the results for the different volume. Fig- 
ure S] and [5] are results of the meson-like correlators and 
their effective masses on a 96 3 x 32 lattice with the same 
quark mass, m q = 0.2, respectively. In the latter section 
(IIIII) . we present simulations in quenched QCD at finite 
temperature, where the finite temperature lattices have 
the similar aspect ratio, N s /N t , and the lowest state en- 
ergy in lattice units to this free quark calculation with 
N s = 96. Although we find the large volume dependence 
between N s — 16 and N s = 96, almost no volume de- 
pendence is seen for N s /N t > 2. At a sufficiently large 
aspect ratio, we find no clear constant contribution, but 
the finite constant contribution exists when the quark 
mass is finite. 



C'(t) 


sinh 






C'(t + 1) 


sinh 





(8) 



Since, however, the overlap of each state is modified by 
a factor of the energy of the state in the analysis, the 
plateau of the effective masses becomes narrow. Further- 
more, in general, a derivative of an observable increases 
the statistical fluctuations, thus the analysis may be diffi- 
cult to apply to actual numerical calculations. Therefore 
we adopt a rather different approach to remove the con- 
stant contribution 2 . We define the midpoint subtracted 
correlator, 

C(t) = C(t)-C(N t /2). (9) 
One can also define the effective mass of the subtracted 



2 We thank Frithjof Karsch for the remark on the midpoint sub- 
tracted correlator. 
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FIG. 6: Effective masses for each channel calculated from 
(left) the midpoint subtracted meson-like correlator and 
(right) the averaged correlator (averaged over the Polyakov 
loop sectors). These are free quark calculations on 16 x 32 
lattices. 



correlator, 



C{t) sinh2 






sinh 2 





(10) 



Figure [6] (left panel) shows the effective mass m s ^ (t), 
defined in Eq. (flQ|) . from the free quark results for 
N s = 16 in Fig. [2l The effective masses are equivalent to 
the usual effective mass shown in Fig. [3] except for the ef- 
fects of the constant mode. The expected energies of the 
lowest two quark states for the S-wave and the P-wave 
states can be calculated from the free quark dispersion 
relation, Eq. ©. The values are shown in Fig. [S] as well 
as in Fig. [3l In contrast to the case of the usual effective 
masses in Fig. [3l the effective masses m |f| (X) approach 
the expected values even in the P-wave states. The anal- 
ysis to avoid the constant contribution works well at least 
in the free quark case. 

The method using the midpoint subtracted correlators 
can also be applied to studies of spectral functions pr{u) 
by a modification of the kernel, 

C{t) = dup r (Lu)K snh (u,t), (11) 
Jo 

K M ~ sinh( W f) ' (12) 

By using the alternative kernel K sn ^°(u),t), it is possible 
to extract the spectral function without the contribution 
in lo <C T from e.g. the Maximum Entropy Method. 

Next we consider the second method to avoid the con- 
stant contribution. The method utilizes the Z 3 transfor- 
mation properties of the correlators, where the Z 3 trans- 
formation is performed on the lattice through the trans- 
formation of temporal link variables, 

U A {x,t) -► e l27m / 3 Ui(x,t), (13) 



where n = 0, 1 and 2, for any x, at a time slice (we 
adopt t = Nf/2). First we consider the Z 3 transforma- 
tion properties for the diagrams in Fig. [Ha) and (b). 
The contribution of Fig. [TJa) has a Z 3 center symmetry, 
while that of Fig. [T^b) is no longer invariant under the 
Z3 transformation. In case of the static limit, one can 
easily understand it by the analogy of Polyakov loop and 
Wilson loop. Even in the non-static case, the properties 
can be derived using a hopping parameter expansion of 
the quark propagator. In any order of the expansion, one 
can show this property like in the static limit. 

Furthermore, for the Z 3 variant contribution, the Z 3 
factor can be factored out. Therefore one can cancel the 
Z 3 variant contribution by averaging the correlators over 
the three sectors which are characterized by the location 
of the Polyakov loop in the complex plane. Here we define 
the averaged correlator C ave (i) from C p0 {t), C pl {t) and 
C p2 (t), which are correlators on the configurations with 
an argument of the Polyakov loop in the complex plane 
fulfilling < arg(P) < 3p ^ < arg(P) < tt and 

— 7T < arg(P) < — respectively, where P is a value of 
the Polyakov loop on the given configuration, 

C aV6 W = \ (C p0 {t) + C pl (t) + C p2 (t)) . (14) 

Each correlator can be calculated on the Z 3 transformed 
gauge configurations. 

The second method to avoid the constant contribu- 
tion is in a way similar to removing the backward quark 
propagations by using a cancellation between results with 
periodic and anti-periodic boundary conditions in the t 
direction. However the second method is beyond a tech- 
nical step to separate the constant mode from the oth- 
ers. The effect of Z 3 symmetrization is to remove from 
the partition function all states with non-integer baryon 
number (whether in quenched or full QCD). This has 
been discussed, e.g., in Ref. [HI • 

Here we should mention that the cancellation of the 
constant mode is not exact in the second method. Be- 
cause, e.g., in the case of a diagram with a 3 times wrap- 
ping quark line, the contribution also has the Z 3 center 
symmetry. Therefore it will not be canceled in the anal- 
ysis. (This situation is also similar to the cancellation 
of the backward quark propagation as mentioned above. 
This cancellation is also not exact. The remaining Z 3 
symmetric contribution can be removed by considering 
the zero-baryon partition function, as discussed as well 
in Ref. [HI-) However, the Z 3 invariant constant con- 
tributions is of 0(e~ 3rrlqLt ) or less, while that for the 
leading (canceled) contribution is 0(e~ mqLt ). Therefore 
the averaged correlator is useful to see a qualitative ef- 
fect of the constant contribution in the correlators, when 
m q L t is not so small. 

From the averaged correlator in Eq. (|14[) . one can de- 
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fine the effective mass m a X- e (i), 



C aYe (t) 


cosh 






c ave o + i) 


cosh 





(15) 



The effective masses m a g- e (t) for the free quark case are 
presented in the right panel of Fig. [6l Similar to the ef- 
fective masses from the midpoint subtracted correlators, 
m a Jj e (£) approach the expected nonzero lowest energies 
even for the P-wave states. However the method does 
not exactly remove the constant contribution, a part of 
the constant contribution remains. 

The method using the midpoint subtracted correlators 
does not work for correlators with finite momentum and 
non-degenerate masses of valence quarks, e.g. heavy-light 
meson. Although these cases are out of the scope of 
this study, the averaged correlator method works even 
for such meson-like operators. Furthermore provided the 
peak structure of spectral function at ui = has a wide 
width to > T due to interactions, the subtracted correla- 
tor method may not even work in the case of correlators 
with zero momentum and degenerate quarks, while the 
averaged correlator method can be applied even in this 
case. 



III. THE CONSTANT CONTRIBUTION IN 
QUENCHED LATTICE QCD AT FINITE 
TEMPERATURE 

A. Lattice setup 

In this section we demonstrate effect of the constant 
contribution and several analyses to avoid it in quenched 
QCD at finite temperature. In the simulation we calcu- 
late the charmonium correlators, which are defined like in 
Sect. UlAl at the charm quark mass. The temporal corre- 
lators (or the spectral functions) of charmonia are impor- 
tant to discuss the J/ip suppression which is considered 
as one of the most important signals of QGP formation 
in heavy ion collision experiments. Several groups have 
investigated the charmonium correlators in finite tem- 
perature lattice QCD 0, H, 0, S Hi ■ In such studies one 
has to take into account the constant contribution to the 
charmonium correlators (or spectral functions) . Because 
this contribution may not be negligible in finite temper- 
ature lattice QCD. In fact it is a genuine effect in finite 
temperature QCD as mentioned in Sect. Ill Al 

In finite temperature lattice QCD the temporal corre- 
lators are restricted within the inverse of the tempera- 
ture. In order to keep sufficient number of data points 
in temporal direction at high temperature, we adopt an 
anisotropic lattice, on which the temporal lattice spacing 
at is smaller than the spatial one a s . The temperature 
is varied by the temporal lattice sizes N t with a fixed 
temporal lattice spacing a t . Such calculations with a 



fixed lattice spacing also have the advantage to investi- 
gate thermal effects in correlators, i.e. one can directly 
compare the correlators at different temperatures with- 
out need for information on the beta-function. In this 
section physical quantities are expressed in at = 1 units 
when there is no comment on the units. 

The gauge configurations have been generated by an 
standard plaquette gauge action with a lattice gauge cou- 
pling constant, f3 = 6.10 and a bare anisotropy parameter 
7g = 3.2108. The definition of the action and parameters 
are the same as that adopted in Ref. [25| . Lattice sizes are 
20 3 x N t where N t = 160 at T = and N t = 32, 26 and 
20 at T > 0. The lattice spacings are l/a s = 2.030(13) 
GeV and 1/at — 8.12(5) GeV which are determined from 
the hadronic radius = 0.5 fm. The physical volume 
size is about (2fm) 3 and temperatures at N t — 32, 26 
and 20 are T = 0.88T C , 1.08T C and 1.40T C respectively. 
The temperatures have been determined by the peak po- 
sition of the Polyakov loop susceptibility from which we 
found the critical temperature around Nt — 28. The fi- 
nite temperature lattices have been generated up to 300 
configurations after 20,000 sweeps for equilibration, each 
configuration is separated by 500 sweeps to reduce the 
autocorrelation. At zero temperature, we have 60 config- 
urations with the same condition. In the calculations, the 
error analysis is performed by the jackknife estimation. 

For the quark fields, we adopt an O(a) improved Wil- 
son quark action with tadpole improved tree level clover 
coefficients. Although the definition of the quark ac- 
tion is the same as in Ref. (25|, we adopt a different 
choice of the Wilson parameter r = 1 to suppress ef- 
fects of lattice artifacts in higher states of charmonium 
[2l| . Then the parameters in the quark action 3 are 
k = 0.10109, c E = 1.911, cb = 3.164 and lF = 4.94. 
Figure [7] shows the usual effective masses calculated with 
the parameters listed above. The results at zero tem- 
perature satisfy a relativistic meson dispersion relation 
with the same anisotropy as the renormalized gauge 
anisotropy a s /a t — 4. The left and right panels in 
Fig. El show the results from the local and spatial ex- 
tended operators. The latter operators are defined by 
Qr(x,t) = J2y 4>(.y)q{x - V,t)Tq(x,t) with a smearing 
function <p(x) in Coulomb gauge. In this calculation 
the spatially extended operators are adopted only in the 
source operators, the sink operators are local in any cases. 
The smearing function is the same as that in Ref. 0], 
i.e. 4>(x) — exp(— ^4|x| p ) where A and P are parameters 
determined by a matching with the charmonium wave 
function as A = 0.2275 and P = 1.258. 

From the zero temperature calculation with the spa- 
tially extended operators we obtain the charmonium 
masses in physical units, m, )c = 3033(19) MeV and 



3 Our lattice setup is very similar to one of Ref. Although our 
notation of the parameters are different from theirs, the corre- 
sponding parameters are consistent with the ones they adopted. 
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FIG. 7: Effective masses from charmonium correlators in 
quenched QCD at zero temperature. The left panel is the 
result with the local charmonium operators and the right one 
is the result with the spatially extended charmonium oper- 
ators. Dotted lines in the left panel show the experimental 
values for each channel. 



FIG. 9: Temperature dependence of the effective masses from 
charmonium correlators in quenched QCD at finite tempera- 
ture. The left panel shows the Pseudoscalar (Ps) channel and 
the right one is the Vector (V) channel. These lowest states 
correspond to r) c and J/ip- 
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FIG. 8: Effective masses from charmonium correlators in 
quenched QCD at finite temperature. These are results with 
local operators for each channel. 

m J/\jj — 3107(19) MeV, which are slightly heavier than 
their experimental values m^f 3 ' = 2980 MeV and 
TO j/>' = 3097 MeV . The experimental values in 
lattice units are also shown in Fig. [71 The results for the 
P-wave states are rather noisy at our current statistics, 
and the actual values are not so important for this study. 
Therefore we do not show the numbers in MeV, how- 
ever their effective masses show a plateau at reasonable 
values. 



B. Finite temperature results 

In this section we discuss the result of the charmonium 
correlators at finite temperature. The definitions for sev- 
eral types of effective masses can be taken over from that 
of the free quark calculations in Sect. [TTJ First we show 
the usual effective masses, m e Q(t) at finite temperature 
in Fig. [H Since these are results obtained with local op- 
erators, there is no plateau region at any temperature or 
channel we chose. It can be expected from zero tempera- 
ture calculation and the effective masses show reasonable 



behavior below T c , at which the charmonium correlators 
does not show large difference from the zero temperature 
results [1, [E II] • We find rather large statistical fluctua- 
tions in the P-wave states, a reason for that is given later. 
Just above T Cl we see drastic changes in the P-wave state 
channels. The changes may cause the dissolution of _v c 
states just above T c as already reported in Ref. 0, 0, HI . 

On the other hand, the S-wave states, i.e. the Ps and 
V channels, show small changes up to 1.4T C . This tem- 
perature dependences are shown in Fig. [9] for each chan- 
nel. When we compare the effective masses at the same 
time slices, we find no large change up to 1.4T C in both 
channels as well. The results are consistent with pre- 
vious lattice studies of charmonium spectral functions 
BSBSUl, in which the spectral functions for the Ps 
and V channels in the deconfinement phase (but at not 
so high temperature) show a peak structure around the 
T] c and J/4> masses like the result at zero temperature. 



C. Midpoint subtracted correlator analysis 

As discussed for free quarks in Sect. HH the constant 
contribution may dominate the meson correlators for 
some channels in the deconfinement phase. In order to 
see the contribution in the meson correlators we calculate 
the effective mass from the midpoint subtracted correla- 
tor, m ^g(t) defined in Eq. (jTUJ) . Figure ITD1 shows the re- 
sults of m^ ) (t) in quenched QCD at finite temperature. 
As we expected, the effective masses at each tempera- 
ture in the S-wave states show small differences from that 
of zero temperature, when comparing them at the same 
time slices. In the P-wave states, however, the drastic 
changes of the usual effective masses m e g(t) are absent 
in the effective masses from the midpoint subtracted cor- 
relators m^g? 3 ^) Furthermore the similar behaviors of 
the effective masses hold till 1AT C as in the case of the S- 
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FIG. 10: Effective masses from the midpoint subtracted char- 
monium correlators in quenched QCD at finite temperature. 
These are results with local operators for each channel. 
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FIG. 11: Time histories of the real and imaginary part of the 
Polyakov loop P at each temperature simulation. 

wave states. The results indicate that the drastic changes 
of the usual correlators (and effective masses) are caused 
only by the constant contribution to the correlators. The 
situation is very similar to the free quark case discussed 
in Sect. HH 



D. Polyakov loop sector dependence of correlators 

Next we discuss the Polyakov loop sector dependence 
of meson correlators as mentioned in Sect. Ill CI First 
of all we show time histories of the Polyakov loop on 
our finite temperature configurations in Fig. 1111 For all 
configurations above T c the Polyakov loop lies in the real 
sector. 

Thus Fig. [5] gives the result on configurations with the 
Polyakov loop being in the real sector. The correlators 
for the imaginary sectors are calculated on configura- 
tions generated by Z% transformations from that in the 
real sector. As an example for the imaginary sectors we 
show the effective masses calculated from C p0 (t) , C pl (t) , 
C p2 (t) and C* ave 0), which are defined in Sect. [EC] in 
FigHH As discussed in Sect. Ill Cl the Z3 transformation 
changes a part of the constant contribution by a factor 
of Z3, i.e. e l2ra7r / 3 where n = 0, 1 and 2. Since our defini- 
tion of the meson correlators is real, the factor acts like 



FIG. 12: Polyakov loop sector dependence of Effective masses 
for the Av channel. pO, pi and p2 sector means the result 
has been obtained on configurations with the Polyakov loop 
restricted to the real and (positive and negative) imaginary 
sectors respectively. Average means averaged over all sectors 
as in Eq. (Tl4)) . 

Re(e i2n7r / 3 ). It means that the constant contributions in 
the imaginary sectors become negative when the value in 
the real sector is positive. Hence the correlators of the 
P-wave states approach a positive constant value like in 
the real sector, as shown in Fig[2l while the correlators 
for the imaginary sectors approach a negative value. As 
a result the effective masses of the P-wave states diverge 
at a time slice as shown in Fig. [T2] 

Here we comment on the large statistical fluctuations 
of correlators for the P-wave states just below T c as can 
be seen in Fig. [8] Below T c , in principle, the constant con- 
tribution vanishes due to the infinite energy of a single 
quark state in the confinement phase. However the con- 
tribution on each configuration may have a small value 
near T c even in confinement phase, like the Polyakov loop 
just below T c in Fig. [Til In this case the effect of the con- 
stant contribution appears in the correlators. Since the 
Polyakov loop sector frequently changes below T c , the 
effects appear as the large fluctuation in the correlators. 

When the averaged correlator defined in Eq. (fT4|) is 
calculated, one can expect that most of the constant 
contributions are canceled like in the free quark case 
(Sect. |HC|) . Figure [lj shows the results of m^ e {t) at 
each temperature. We find almost the same behavior 
as in Fig. [10] The results also suggest that the drastic 
changes in the P-wave states just above T c are coming 
from the constant contribution only. 



E. Results with spatially extended operators 

Finally we present results of meson correlators with 
spatially extended operators. The spatially extended op- 
erators are constructed with a smearing function which 
yields a spatial distribution of quark (and anti-quark) 
source (s). The distribution also changes the momentum 
distribution of the quark propagation. For example, a 
point source function yields any (lattice) momenta with 
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FIG. 13: Effective masses from the averaged charmonium cor- 
relators (averaged over the three sectors of the Polyakov loop) 
in quenched QCD at finite temperature. Configurations with 
the Polyakov loop in the imaginary sectors have been obtained 
by a Zz transformation of the corresponding real sector con- 
figuration. These are results with local operators for each 
channel. 



FIG. 15: The same as in Fig. [10] but for spatially extended 
operators in each channel. 
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FIG. 14: The same as in Fig. [8] but for spatially extended 
operators in each channel. 



equal distribution, while a wall source function yields 
only quark propagation with zero momentum. In the 
case of our source the (smearing) function provides quark 
propagation with any momenta but quark propagators 
with small momentum are enhanced. Therefore the cor- 
relators with spatially extended operators defined, as in 
Sect. IIII A[ yield (relatively) large overlap with the lower 
energy states. 

Figure [JJ] is the result for the usual effective masses 
with the spatially extended operators. Other conditions 
are the same as in Fig. [5] Below T c , in contrast to the case 
of local operators the effective masses reach a plateau 
due to the larger overlap with the lowest state. Above 
T c the effective masses of the P-wave states change more 
than in the case of local operators. Because the constant 
contribution is enhanced by the smearing function more 
than the other contributions in the case. The effect of 
the constant contribution is also present in the V chan- 
nel. At least in the free quark case, as we discussed in 
Sect. Ill A[ the constant contribution to the V channel is 
smaller than that to the P-wave states for finite quark 
masses and small volumes. In fact, even in quenched 
QCD, it is difficult to see the effect in the V channel in 
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FIG. 16: The same as in Fig. [13] but for spatially extended 
operators in each channel. 



local operators. However we find small but finite changes 
even in the V channel above T c while no changes are seen 
in the Ps channel as we expected. 

We can say that the small changes in the V channel is 
caused by the constant contribution, because the analysis 
with the midpoint subtracted correlators shows almost 
no changes in the V channel above T c as it can be seen 
in Fig. [121 Similar to the results with local operators, 
the effective masses of the P-wave states also show small 
change at any temperatures. 

Figure [TBI is the same figure as Fig. [Jj2 but with spa- 
tially extended operators. We find similar behavior for 
both cases of operators. 



IV. DISCUSSION 

Let us discuss details on \c states using the results 
of the previous sections. In this study we find that the 
drastic changes of the (usual) correlators in the P-wave 
states are caused by the constant contribution, while the 
correlators without the constant mode yield small change 
till, at least, 1.4T C even in the P-wave states. The changes 
are of similar size as that of the S-wave states as one can 
see in Fig. [17] and Fig. [18] 

In the case of the spatially extended operators, the av- 
eraged effective mass m^g- c (t) for the Av channel grad- 
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FIG. 17: Temperature dependences of the effective masses 
from the midpoint subtracted charmonium correlators (left) 
and the averaged charmonium correlators (averaged over the 
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at finite temperature. These are results with local operators 
of the V and the Av channels, whose lowest states correspond 
to the J/ip and Xci states. 
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FIG. 18: The same as in Fig. [17] but for spatially extended 
operators in the V and the Av channels. 



ually reduces as the temperature increases. The change 
seems to be effected by the constant contributions with 
the Z 3 center symmetry as discussed in Sect. Ill Cl In fact, 
we can not find any hints of such a change in m 8 ^ 3 ^) in 
the left panel of Fig. [T5I Furthermore it is reasonable to 
find the effect only in case of the spatially extended op- 
erators, because the Z 3 symmetric constant contribution 
is enhanced by the smearing function. 

Although the constant mode coming from the scatter- 
ing process is important to investigate some transport 
phenomena, it does not affect the spectral function in 
the uj 3> T region for the charmonium correlators with 
zero spatial momentum. Therefore the dissolution of y_c 
states just above T c as discussed in e.g. Ref. 0, 0, Q 
might be misleading. When there is the constant contri- 
bution, the hadronic state is no longer the lowest state. 
The situation provides some difficulties in the analysis 
to reconstruct the spectral function in the higher energy 
part uj > T even if the MEM is applied [27j]. The difficul- 
ties may be understood by the fact that most studies of 
the charmonium spectral function can not reproduce the 
excited states correctly especially at finite temperature 



(even below T c ). Furthermore, in general, correlators for 
the P-wave states are noisier than that for the S-wave 
states as one can see in Fig. 

In principle the MEM analysis can extract the correct 
spectral function even if the constant contribution exists, 
however it is difficult to reproduce the non-lowest part of 
the spectral function correctly by the MEM analysis at 
finite temperature with typical statistics. In order to in- 
vestigate the properties of Xc states above T c , a careful 
analysis taking the constant contribution into account is 
necessary. A correct MEM analysis of the correlator with 
and without the constant contribution should give the 
same results, up to the low energy region uj <C T. Such 
a comparison can be easily performed by the MEM anal- 
ysis with the alternative kernel of Eq. (fT3i)) for midpoint 
subtracted correlators. Failure to do so signals unreliable 
MEM results. 

We also find a Polyakov loop sector dependence of the 
constant contribution in quenched QCD. To study trans- 
port coefficients from the behavior of spectral function 
at uj — one has to pay attention to the Polyakov loop 
sector especially in quenched QCD calculations. From a 
practical viewpoint, the dependence on the Polyakov loop 
sector is useful to estimate the constant contribution as 
discussed in Sect. Ill CI and Sect. lIIIDl 

The constant contribution may also explain some other 
results of previous studies of temporal correlators. First, 
in the study of the spectral function at finite tempera- 
ture by the MEM analysis, one sometimes find a kind of 
divergence of the spectral function near uj — 0, it may be 
caused by the constant contribution. Of course it may 
be a simple failure of the analysis, however the contribu- 
tion near ui = ought to exist except for the Ps chan- 
nel. Secondly, in the S-wave states of charmonium, sev- 
eral groups have reported that the correlator for the Ps 
channel shows almost no change above the deconfinement 
transition, while that of the V channel shows a small but 
significant change [1, H, [28| . This can also be seen from 
our calculation in Sect. IIIII It can be explained by the 
constant contribution as we did in this paper. Thirdly, 
when one considers a spatial boundary condition depen- 
dence of the charmonium spectrum at finite temperature 
as in Ref. [28| , the P-wave states (in the free quark case) 
have the lowest energy using a "hybrid boundary" , which 
is an anti-periodic boundary in one spatial direction and 
periodic in the other directions. On the other hand, the 
P-wave state has a larger energy with periodic boundaries 
in any spatial directions. However the (usual) effective 
mass of the P-wave state correlators for free quark cal- 
culations provides a smaller value in periodic boundary 
conditions than in the hybrid boundary case. This un- 
expected behavior of the effective mass is also explained 
by the scattering processes because of suppression of the 
scattering contribution by the anti-periodic boundary in 
a spatial direction. Finally, the Polyakov loop sector de- 
pendence of the chiral order parameter in quenched QCD 
|2g | may be also explained by the effect of (Z 3 variant) 
wrapped contributions in the chiral condensate. 
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As we discussed in the last section, the constant con- 
tribution affects the temporal correlators at finite tem- 
perature in various cases. Therefore it is important to 
take this contribution into account for studies of tempo- 
ral correlators in lattice QCD at finite temperature. 
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